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Abstract 

Due to the smallness of the electromagnetic coupling, photons escape from the hot and dense 
matter created in an heavy-ion collision at all times, in contrast to hadrons which are predominantly 
emitted in the final freeze-out phase of the evolving system. Thus, the thermal photon yield carries 
an imprint from the early evolution. We suggest how this fact can be used to gain information 
about where between the two limiting cases of Bjorken (boost-invariant expansion) and Landau 
(complete initial stopping and re-expansion) hydrodynamics the actual evolution can be found. We 
argue that both the rapidity dependence of the photon yield and photonic HBT radii are capable 
of answering this question. 

PACS numbers: 25.75.-q 
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I. INTRODUCTION 


Q 


The boost-invariant hydrodynamic model proposed by Bjprken [l| for the description of 
ultrarelativistic heavy-ion collisions is frequently used at RHIC energies for estimates of the 
initial energy density in heavy-ion collisions or the lifetime from the measured Hanbury- 
Brown Twiss (HBT) correlation radius Ri on g Q as well as in hydrodynamical descriptions 


of the evolving system (see. e.g. |3J). 

While the original notion of boost-invariance is an asymptotic concept, its application to 
RHIC energies usually implies two things: 1) the distribution of matter in some finite interval 
around midrapidity is assumed to be (almost) independent of rapidity and 2) the longitudinal 
dynamics is assumed to be unaccelerated expansion which in turn means that momentum 
rapidity y = | In is always equal spacetime rapidity r) s = \ In (this is not true in 
the presence of longitudinal acceleration). 

n contrast, charged meson rapidity distributions as obtained by the BRAHMS collaboration 




do not show a flat plateau around midrapidity even at top RHIC energy. The distributions 


are however well described by Landau hydrodynamics |5| as argued in [6|. Likewise there 
is no boost invariance seen in the rapidity dependence of elliptic flow as measured by the 
PHOBOS collaboration [7]. 

In a model framework adjusted to reproduce the full set of observables characterizing the 
hadronic freeze-out, i.e. single particle transverse mass spectra and rapidity distributions 
and two particle HBT correlation radii jg| it was found that simultaneous agreement with all 
data sets can only be achieved if the assumption of a boost-invariant expansion is dropped. 
In fact, a sizeable difference of Ay = 2-1.8 between initial and final width of the source in 
momentum space rapidity is required. 

This, however, is rather indirect evidence since it rests on a backward extrapolation of the 
observed final state. In contrast, thermal photons would offer the opportunity to test the 


longitudinal evolution directly 9]. The essential idea is as follows: In a Landau scenario, 
the source is initially very narrow around midrapidity. Since the hard photon emission rate 
is strongly temperature dependent, the dominant contribution to the photon yield arises 
from early times. Thus, we expect that the hard photon yield as a function of rapidity 
shows a thermal smearing of the initial (narrow) source extension in rapidity. In contrast, 
in a boost-invariant expansion we expect a much broader distribution reflecting the initial 
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distribution of matter across a large rapidity interval. 

There is an additional factor which needs to be taken into account: Due to its large initial 
extension in rapidity, a Bjprken scenario leads to much more rapid cooling than a Landau 
one. Hence, while in a Landau scenario the hot, early phase will be dominant, this is not so 
in a Bjprken framework. The different weights of the contributions of early times and late 
times are expected to leave a characteristic imprint on HBT correlation radii measured even 
at midrapidity. 

In this work, we discuss both ideas and demonstrate what predictions for the photonic 
observables can be made using either the scenario determined from a fit to spectra and HBT 


m 


or a Bjprken or a Landau one. 


II. THE MODEL FRAMEWORK 


Several calculations studying photon emission based on a hydrodynamical fireball evolution 
model have been made so far for different collision systems and energies jltj, [ill . fli! [lij . In the 
present study investigating 200 AGeV AuAu collisions, we will use a parametrized evolution 
model instead which allows for a complete description of hadronic transverse mass spectra as 
well as HBT correlation parameters 8] and which can easily be tuned to interpolate between 
Bjprken and Landau dynamics. 

The model for the evolution of hot matter is described in detail in 
present the essential outline and focus on (almost) central collisions: 

For the entropy density at a given proper time we make the ansatz 


14j . Here we only 


s(r, r] s , r ) = NR(r, r) • H(rj s , r) 


( 1 ) 


with r the proper time measured in a frame co-moving with a given volume element and 
R(r, r), H(r ] s , r) two functions describing the shape of the distribution and N a normalization 
factor. We use Woods-Saxon distributions 


R(r,r ) = 

H(v s , t) 


= 1 / 


1 + exp 
I 1 + exp 


Rc{r) 


d v 


Vs ~ H c (t) 
Vws 


( 2 ) 


for the shapes. Thus, the ingredients of the model are the skin thickness parameters d ws 
and i] ws and the parametrizations of the expansion of the spatial extensions R c {t), H c (t ) as 
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a function of proper time. From the distribution of entropy density, the thermodynamics 
can be inferred via the EoS and particle emission is then calculated using the Cooper-Frye 
formula. For simplicity, we assume that the flow is built up by a constant acceleration a±, 
hence R c (t) = R° c + yyT 2 with R° c an initial radial extension as found in overlap calculations. 
The rapidity distribution is assumed to grow from some initial width 2 ■ x/q to a final width 


2 • yp. This determines the extension of the emitting source in spacetime rapidity r/ s 


8 , 


H]. 


In |8], the model parameters have been adjusted such that the model gives a good description 
of the data. This implies an initial rapidity width of y 0 = 1.7. In order to compute a Bjprken 
scenario, we set the initial width of the rapidity distribution equal to the final distribution 
width T/o = y F . For a Landau scenario we choose y 0 = 0. In both cases we readjust the 
model parameters such that the single particle spectra are reproduced (this implies losing 
agreement with the HBT data). 

The spectrum of emitted photons can be found by folding the photon emission rates for 
the quark-gluon plasma (QGP) phase [3 and for a hot hadronic gas Q with the fireball 
evolution. In order to account for flow, the energy of a photon emitted with momentum 
k M = (kt, k t , 0) has to be evaluated in the local rest frame of matter, giving rise to a product 
with r) the local flow profile. Following the results in we assume for the 

spatial dependence of the flow field the relations y = f{j) ■ q s and y± = g(r) ■ r with y L the 
transverse rapidity and /, g two functions determined by the evolution. The distribution 
of entropy density is manifest in the dependence of the temperature T = T(g s ,r,r ) on the 
spacetime position. In order to account for the breakup of the system once a temperature 
Tp is reached, a factor 9{T — Tp) has to be included into the folding integral. 

Using the folding integral of the rate with the fireball evolution as emission function S(x, K ) 
(describing the amount of photons with momentum emitted at spacetime point x we 
calculate the HBT parameters as [3,13] 


R lde = ( V 2 ) R lnt = {(2 ~ P±tf) R \ong = (?) 


(3) 


with = Xf, 


{xp) and 


{f(x))(K) = 


f d 4 xf(x)S(x,K) 
f cl 4 xS(x, K ) 


(4) 
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III. RAPIDITY DEPENDENCE OF HARD THERMAL PHOTON EMISSION 


We show the resulting spectra of hard thermal photons in the momentum range between 1 
and 4 GeV in Fig. |T| at two different rapidities. 




FIG. 1: The hard thermal photon spectrum at midrapidity 
(y=2, right panel) for the best fit scenario described in 8], 


(y=0. left panel) and forward rapidity 
a Bjprken and a Landau scenario. 


It is instructive to observe that both slope and absolute yield changes strongly as a function 
of y for the Landau scenario. This reflects the fact that the initial high temperature phase 
(leading to a relatively flat slope) never radiates out into the y = 2 slice — only in the 
later stages when hot matter expands across y = 2 there is a significant contribution, albeit 
from matter with a much lower temperature, leading to a steeper spectral slope and reduced 
yield. 

In contrast, the photon yield from a Bjprken scenario is practically unchanged as a function 
of rapidity, reflecting the approximate boost-invariance. 

In order to highlight the differences more clearly we show in Fig. [21 the fc^-integrated yield (1 
GeV < kx < 4 GeV) at rapidity y 0 divided by the integrated yield at midrapidity. This choice 
has the additional advantage that model-dependences such as the precise normalization of 
the emission rates tend to cancel out. 

The different longitudinal source structure is now directly apparent. The Landau scenario 
is characterized by thermal smearing of about 1 unit of rapidity of a source at midrapidity 
(without any longitudinal flow) whereas the Bjprken scenario shows the broad distribution 
of matter across ~ 3 units of rapidity at all times. A measurement of the thermal photon 
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FIG. 2: Integrated yield (1.0 GeV < kx < 4 GeV) of thermal photons as a function of rapidity j/o- 

yield at midrapidity and at y 0 — 2 would well be capable of making a distinction between 
the three scenarios. 


IV. HARD THERMAL PHOTON HBT AT MIDRAPIDITY 


HBT correlation measurements do not measure the true geometrical size of the source but 
rather a region of homogeneity 17, Q] which is only identical with the geometry for vanishing 
flow gradients in the source. For finite flow gradients, the measured correlation radii show 
a characteristic falloff with the correlated pair momentum k x- The precise shape of the 
correlation radii as a function of transverse momentum results from a complex interplay 
between temperature and flow during the whole evolution. 

Nevertheless, we can formulate some basic expectations. Due to the high initial compression, 
the hard photon yield from a Landau scenario is expected to be dominated by the initial 
phase of the expansion. In this phase, however, there is no significant transverse flow (which 
builds up gradually driven by transverse pressure) and the geometrical size of the source 
in longitudinal direction is very small (for complete stopping it is given by the Lorentz- 
contracted size of the overlapping nuclei). Thus, we would expect only a weak falloff of R S id e 
with kx and Ri on g to be determined primarily by the spatial resolution scale of photons with 
a given momentum. 

In contrast, a Bjprken expansion may well receive significant relative contributions to the 
yield from later stages due to the shorter duration of the inital hot phase. This would imply 
a slightly larger R S id e for vanishing kx but a stronger falloff with kx and an increased value 
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Riong as compared to the initial size at equilibration time. The relevant underlying scale for 
Rside is in all cases the nuclear overlap radius. 




FIG. 3: Hard photon HBT correlation radii R S ide (left panel) and Ri on g (right panel) for the best 
fit scenario described in [s | , a Bjprken and a Landau scenario. 

The result of the calculation can be seen in Fig. El To a good degree, the expected behaviour 
is indeed seen. In particular, the different falloff of R S ide for kx > 2.5 GeV appears to be a 
good indicator of the longitudinal dynamics. Ri ong in contrast is presumably only capable 
of identifying a scenario very similar to a Landau one, otherwise the qualititive behaviour 
of the different curves is too similar. Note that the observed Riong for the Landau scenario 
could not be as small as shown in the plot due to constraints posed by the uncertainty 
relation which doesn’t allow to narrow down the photon emission region to arbitrary small 
size. 


V. THE ROLE OF PRE-EQUILIBRIUM PHOTONS 


It is well known that in addition to thermal photons a sizeable contribution of prompt 


lioton yield, 
d,M). This 



photons (calculable in perturbative QCD) is expected to contribute to the hare 
and various attempts have been made to calculate its magnitude (see e.g. [Hj, 
contribution might well outshine the signals proposed here and change the conclusions. 

In order to address this question carefully, we do not only have to take into account the 
primary hard scattering processes as a potential source of photons but also hard re-scattering 
processes as the system approaches equilibrium. Therefore we use here the VNI/BMS parton 
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cascade model (PCM) to estimate the role of pre-equilibrium hard photon production 
23]. 
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There is still the caveat that the re-scattering described by the PCM does not lead to a 
Landau-like stopping of the incoming matter, nevertheless we use the results to gain some 
intuition in the orders of magnitude involved. 

Including the LPM suppression in the PCM, we find that thermal photons may dominate 
the yield below 2-2.5 GeV for the best fit and the Bjprken scenario whereas they would 
dominate the yield in the whole momentum range for a pure Landau evolution 24]. 

Since the photon yield drops (almost) exponentially with kt, this implies that the rapidity 
dependence of the integrated yield would still be a reliable signal (being dominated by the 
low k T yield). 

However, the behaviour of the HBT correlation radii in the interesting region above 2 GeV 
is likely to be distorted by pre-equilibrium photons (which would incidentially resemble 
Landau dynamics as they are characterized by small transverse flow). 

Turning the argument around, a simultaneous measurement of the HBT correlations at 
midrapidity and of the integrated yield at forward rapidity could still provide valuable 
insight into the magnitude of the pre-equilibrium contribution at different momenta. A 
detailed investigation of this question is however beyond the scope of this work. 


VI. SUMMARY 

We have argued that photons provide a direct measurement of the early longitudinal dy¬ 
namics of a heavy-ion collision which can otherwise only be inferred indirectly from hadronic 
probes. The underlying reason is that due to the smallness of the electromagnetic coupling 
the measured photon yield represents an integral over the whole fireball evolution rather 
than a snapshot at breakup. 

In particular, we have argued that the rapidity dependence of the hard photon yield is a good 
probe to distinguish between Landau and Bjprken-likc dynamics since it directly reveals the 
rapidity extension of the emission source. Since we compare the rapidity dependence of 
a ratio of integrated yields many uncertainties associated with the calculation of emission 
rates drop out and the result mainly reflects kinematic properties of the source. 

In addition, we have investigated the potential of using HBT correlation measurements at 





midrapidity to investigate the longitudinal evolution. HBT correlations show what part 
of the evolution dominates the photon yield rather than directly reflecting longitudinal 
dynamics. We found that the falloff of R S id e with k t above 2.5 GeV would indeed give a 
good indication if the photon emission is dominated by matter without significant transverse 
flow or not, however this signal is easily obscured by pre-equilibrium photon emission which 
would never show significant transverse flow. 

Nevertheless, both measuring the rapidity dependence of the hard photon yield and the 
HBT correlation parameters at midrapidity are capable of revealing interesting details of 
the early fireball evolution which cannot easily be studied otherwise. 
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